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We calculate the dilepton production rate from a thermalized quark-gluon 
plasma in heavy-ion collisions at RHIC energies. Higher-order QCD correc- 



tions are included by using an analytical correction factor K®, which takes 



into account the qq initial-state interactions. We show that the analytic cor- 
rection factor gives very good agreement with experimental Drell-Yan data 
and leads to large enhancement of the thermal dilepton emission rates. We 
compare the thermal dilepton yields with the expected production from open- 
charm decays and Drell-Yan background and assess the prospects of observing 
thermal dileptons from the quark-gluon-plasma at invariant masses of a few 
GeV. 
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High-energy heavy-ion collisions have become the focus of intense experimental and the- 
oretical research in recent years because of the possibility of producing hadron matter in 
a deconfmed quark-gluon plasma (QGP) state during such collisions ]!]]. Various signals 
have been proposed as a probe of this new phase of nuclear matter. In particular, electro- 
magnetic signals comprising dileptons and photons, are recognized as direct probes for the 
QGP phase 0-fl], since they weakly interact with the hadronic medium in which they are 
produced. The magnitudes of the dilepton yields depend on the cross sections for the basic 
process producing the dileptons. 

The dilepton signals indicative of the QGP phase are those produced by the annihilation 
of qq pairs through virtual photon intermediate states. To be meaningful signals, they must 
be clearly delineated from other sources of dilepton pairs like Drell-Yan process, open charm 
decays, vector meson decays, and so on. The problems associated with these backgrounds 
have been discussed in the literature |5|-|7|,£|,|8||| and their contributions are known to vary 
with the invariant masses of the dileptons and the kinematic regions under consideration 
PJiTi|-|r5[ . For dilepton pairs in the invariant mass region of 2 to 7 Gev, the charm dileptons 



and DY dileptons are expected to be the major background to detection of thermal dileptons 
from the quark-gluon-plasma phase. We shall discuss our thermal dilepton predictions in 
the light of these backgrounds. 

Thermal dileptons as well as Drell Yan pairs originate from the electromagnetic annihila- 
tion of quark-antiquark pairs through intermediate virtual photons. While the surroundings 
that host the processes for dilepton production from the thermal plasma or the DY pro- 
duction during nuclear collisions are vastly different, the basic process corresponding to the 
channel qq — > 7* — > is the same for the two cases. For comparative analyses, there- 
fore, the dileptons should be corrected similarly for the effects that modify the tree-level 
diagrams, whether the dilepton source is the thermal plasma or the DY mechanism. For 
the latter processes, it has been already generally recognized in perturbative QCD that the 
lowest-order Feynman diagrams give only an approximate description [1^-^]. In order to 



bring the lowest-order QCD predictions into agreement with experiment, the lowest-order 



DY results must be multiplied by a phenomenological i^-factor, which has a magnitude of 
the order of 3. The i^-factor can be accounted for in terms of higher-order QCD corrections, 
the most important effects being the vertex correction arising from the initial-state interac- 
tions of the qq, prior to their annihilation into the virtual photon producing the dileptons 
|n] , |r5|] . These higher-order QCD corrections are operative also for qq annihilation in the 
QGP, and the dilepton production cross section in the QGP must be similarly corrected. 

In this work we investigate dilepton production from the QGP in nucleus-nucleus colli- 
sions at RHIC energies, by incorporating these QCD corrections through an analytic cor- 
rection factor presented recently ||19|| . We compare our calculated rates with those expected 
from open charm production and the DY process in the energy range of a few GeV. Our 
calculated values of open charm dileptons exceed those of the DY process, in conformity 



with recent findings [|12j. However, we show that the thermal dilepton yield, computed 
with the QCD corrections can be comparable and in some scenarios can be well above the 
background, so that prospects for its observation at RHIC energies becomes promising. 
We use the analytical correction factor for the basic reaction processes involving q and 



q pi, for the whole range of relative energies. For low energies, the correction factor is 
obtained by studying the distortion of the wave function by virtue of the q-q color potential 
arising from virtual gluon exchange. At high energies, the correction factor is chosen to 
match well-known PQCD results which give good agreement with the experimental correc- 
tion factor for the initial- or final-state interactions as applicable. An interpolation to join 
the correction factor from low energies to these well-known perturbative QCD correction 



factors at high energies is made by following the procedure suggested by Schwinger f20 |. 

The analytical correction i^-factor for the annihilation or the production of a qq pair 
with an invariant mass M 2 = s is explicitly given by [[19 



where the flavor label q in K^ l \q) is included to indicate that K^i depends on the quark 
mass m q , the superscripts (i) and (/) denote qq initial-state annihilation and qq final-state 



production respectively. 
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In the above equations, v is the relative asymptotic velocity for the quark and the antiquark 
in their center-of-mass system 



^2_ 4sm 2)l/2 

s — 2m q 
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and a e g is the effective strong interaction coupling constant related to the strong interaction 
coupling constant a s by the color factor Cf 



«eff = C f a s . 



For dilepton production, the q and q must be in the color singlet state, i.e. Cf 



(5) 



■|, and 



the running coupling constant is taken as |T6|,|T7 
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To demonstrate the validity of the correction factor, we apply it first to the Drell-Yan 
case and compare with known experimental data. 

The lowest-order Drell-Yan distribution for nucleon-nucleon collisions is given by 

„2 ,„ \ 2 



dMdy ~ 3sN c M , 

a L q=u,d 



J2 (%) [q a q (xe»)q b q (xe~y) + q a q (xey)q b q (xe~y)} 
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where x = M/yfs and y is the rapidity. q q ,b (xe y ) and q q ,b (xe y ) refer to the quark and 
antiquark distributions in the nucleons a and b respectively. 
Correcting for initial-state color interaction we obtain 



dMdy 3sN c M ^ 



E ^%)(-) [q a M y )qt(xe-v) + F q (xey)q b Jxe-v)} 
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The DY distributions for equal-nuclei nucleus-nucleus collisions are obtained from the 
nucleon-nucleon case by using 



rt 2 N AA 3 d 2 rr NN 

u Jv fly _ ° ^4/3 " DY /q\ 
dMdy ~ An(r' ) 2 dMdy ' U 

where A is the atomic number of the colliding nuclei and r' = 1.2 fm. 

Figure [l| shows the FNAL-605 experimental data plj for the differential cross section 

sd 2 a / d^/rdy as a function of y/r = M/^fs for different rapidity y intervals, compared to 

the corrected cross-sections using the Duke and Owens structure functions with A = 0.2 



GeV ||22|| ) obtained by multiplying the lowest-order Drell-Yan calculations by the correction 



factor if". No further multiplicative factors are needed to obtain the very good fit shown in 
fig. 0, which thus demonstrates the reliability of the correction factor It is interesting 
to note that the function from Eq. (||) which gives a good fit in Fig. [I] is not a constant 
of the invariant mass M. It is equal to 2.5 for M = 3.9 GeV and 1.8 for M = 17 GeV. There 
is thus a 30 % variation of K^ 1 ' as M varies from 4 to 17 GeV. 

For lowest-order calculations, the rate for the production of dileptons with an invariant 
mass M per unit four-volume in a thermalized quark-gluon plasma (with three flavors) 
depends on the temperature T of the system and can be written in the form ]7|,[l| 



<"W ~ NcN} E < A - '4 ™* (%) , (!0, 



dM 2 d 4 x s ^ \ej 2(2tt) 4 V M 2 L \T 
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where <J q (M) is the lowest-order qq — > cross section at the center-of-mass energy M 
given by fll 
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mj and m q are the rest masses of the lepton / and the quark q respectively, and K\ is the 
modified Bessel function of first order. 

The integration of Eq. ( |T(]| ) over the transverse dimensions and the proper time interval 
ro (the formation time) to r c (the hadronization time) gives the invariant mass-squared 
distribution of dileptons produced from the system while it was in the quark-gluon plasma 
phase as |7],[lJ 

/ » \ 2 ~ { A/r\ / /lm 2 \ 7 Q^2rpQ r / ]\/T\ / 1\/T\ 1 
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where A is the transverse area of the collision region, To and T c are the initial and final 
temperatures of the ideal QGP phase, and 

H{z) = z 2 (8 + z 2 )K {z) + 4z(4 + z 2 )K x {z) . 

To delineate the QGP dileptons from other sources of dileptons, particularly the Drell- 
Yan background, we need to have good estimates of the dilepton production from the dif- 
ferent sources at matching levels of accuracy. It is hence necessary to modify Eq. (P) by 
incorporating the higher-order QCD corrections. The interpolation from low to high relative 
velocities extends the usefulness of the correction factors given by Eq. ([!]) to account 
for higher-order QCD effects. The modified rate of dilepton production then becomes 



dMHy - M ^^HT) & M*J 



■ (13) 



T J \T C 

To calculate Eq. (|12|) or Eq. (|i~3|), we need to know the temperatures T and T c , and 
the "formation" or "materialization" time To. We use T c = 180 MeV. The Monte Carlo 
simulation program "MARCO" |]23[| , which is based on the Glauber multiple-collision model 
and reproduces the peak value of dN/dy in nucleus-nucleus collisions, is used to calculate 
the rapidity density of the produced particles (mostly pions). The peak energy density e is 
related to the peak rapidity density dN/dy\ pea k in the Bjorken scaling hydrodynamic model 



m T dN 



t A dy 



(14) 

peak 



where m T = ^jm\ + p^, is the pion mass, and p T is the pion transverse momentum, 
which varies with the center-of-mass energy of the colliding nuclei as ||T| 

p T = 0.27 + 0.037 ln(Vs). 



The temperature T is then approximated as P^jJUl 



30 e 
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The value of the formation time To is found to play an important role in determining the 
value of the initial temperature T of the QGP. It has been given different values in the 
literature. In the following, we use the values obtained from MARCO for the peak dN/dy, 
typically 2144 and 231, in the central collisions of Au + Au and S + S, respectively, at RHIC 
energies. The common argument assumes that r should be of the order of a typical strong 
interaction time scale of 1.0 fm/c, as proposed by Bjorken [pi]] . Using this value, Eq. (|i~4"l) 



and Eq. (|T5D then give T = 262 MeV for Au+Au and T = 202 MeV for S+S. Shuryak et 
al. H25l , |26l assume a different scenario for calculating the thermalization time tq. In the "hot 
glue" scenario, the thermalization time of the quark-gluon plasma is expected to be within 
the range < r < 0.3 due to gluon rescattering and production mechanisms. This picture 
gives a temperature much higher than that of the Bjorken scenario. Ruuskanen 0] assumes 
that the quark-gluon plasma may thermalize at an early time r = 0.5 fm/c which gives 
T = 312 MeV for Au+Au and T = 241 MeV for S+S. Alternatively, Kapusta uses a 
different approach to estimate To. Assuming that the partons will have a thermal distribution 
right after they are put on the mass shell, and that the initial rapidity density is the same as 
the observed one, Kapusta concludes that r T = constant. Using the uncertainty principle 
and putting AE = 3T , he then estimates the constant to be 1/3, i.e 

n>T = l/3, (16) 

which gives T = 590 MeV and r = 0.11 fm/c for the collision Au+Au and T = 421 MeV 
and r = 0.16 fm/c for the collision S+S. 

Due to the importance of the charm dilepton background, it is necessary to estimate 
the charm contribution for a comparitive study. For charm production, the gluon fusion 



channels as well as the qq annihilation via virtual gluon modes must be included ||27|| . 

The basic cross section for qq — > g* — > cc, averaged over initial and summed over final 
colors and spins can be written as 



where rj = 4m^/M^ g ., with m c being the mass of the charm quark and M c5 being the invariant 
mass of the produced cc pair. 

The corresponding expression for the gluon fusion mode, averaged over initial gluon types 
and polarizations and summed over final colors and spins is [I] 

'-«■> ^ + ' + K^) - (I + 1")^) } 

The above cross sections must be convoluted over the quark and gluon distributions respec- 
tively as for the DY case to obtain the overall charm yield for comparison with the QGP 
signals. The charm production can then be written for the nucleon-nucleon case as 

j2„NN 



d o c - c 



dMc _ cdyc _ c K[ ( £ a q ,(M c5 )[q a q (xey)q b q (xe~y) + g^xe^V^f 

+ ([g a (xey)g b (xe-y)} + [a «-> bty m {M*)} (19) 
where x = M c5 /^, M c5 , y c5 are the invariant mass and the rapidity of the cc pair. The 



factor K in this case is an empirical factor and is usually taken to be 3 |28 . 

The cc pairs produced by the qq or gg reactions hadronize to D mesons which subse- 
quently decay through semi-leptonic channels, such as D — ► I + X to leptons. Leptons from 
a decaying DD pair comprise a dilepton and the charm production rates can be mapped 
into the resulting decay dilepton distributions with invariant mass Mi + i_ and rapidity yi + i- 
by a Monte Carlo technique to generate the dileptons from charm decays. 

The nuclear-nuclear charm dilepton cross section can be obtained from the nucleon- 
nucleon one according to 

d 2 Nj+f- 3 ^ 4/3 d 2 af+f- 



dM l+l -dy l+l - 47r(r^) 2 dM l+l -dy l+l - 
The dilepton production rates per unit time from a thermalized QGP, which may be 
formed after the collisions of Au+Au and S+S at an energy = 200 GeV, for three 
different values of the variable To are shown in Fig. |2] and Fig. |^, respectively. The solid 
curves in Fig. g(a) and Fig. |3|(a) show those rates given by Eq. ( TT2] ) using the results of 
the lowest-order cross section of dileptons production. The solid curves in Fig. ||](b) and 
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Fig. |](b) show the same rates taking into account the correction factor K^> [Eq. ([13])]. The 
dotted-dash curves in both Fig. |2| and Fig. |3| show the Drell-Yan rates. The dotted curve is 
the dilepton rates from the open-charm mesons decays. 

The virtual photon annihilation mode selects the annihilating qq pair to be in the color- 
singlet states and thus the interaction between q and q, before they annihilate, is attractive. 
The effect of qq initial-state interaction leads to an enhancement factor as shown in 
Fig. |2] and Fig. Because the coupling constant increases and the relative velocity decreases 
as the invariant mass decreases, the correction factor rises considerably with the decrease of 
invariant mass. The effect of the higher-order QCD corrections is to enhance the tree-level 
dilepton cross section by a factor of about 5 at y/s = 1 GeV, and by a factor of about 3 at 

= 2 - 3 GeV. 

It is seen that for the collision Au+Au, the temperatures T = 262 MeV and T = 312 
MeV lead to the result that the dilepton rates from the Drell-Yan process (the dashed curves 
in Fig. 0) exceed those from the QGP (the solid curves in Fig. 0(a)) above 1 GeV, and the 
charm dilepton contribution is even higher. The temperature To = 590 Mev allows the 
uncorrected thermal rate to exceed the DY and charm dileptons up to almost 6 Gev and 
indicates it may be possibile to detect the signal. However, the inclusion of the correction 
factor changes the situation dramatically. It increases the thermal rates (solid curves in 
Fig. 0(b)) so that they exceed the Drell-Yan background up to an invariant mass of M = 2 
GeV even for the lower temperature models and is comparable to the charm dileptons, which 
may prove to be detectable if the contribution of the charm dileptons can be found by inde- 
pendent charm measurement. For the temperature To = 590 MeV, it further improves the 
prospects to detect the QGP dileptons since their rates are now well above the background 
up to M = 7.5 GeV. 

Similarly, for the collision S+S at an energy y^i = 200 GeV, as shown in Fig. ^| at lower 
temperatures T — 202 MeV and T — 241 MeV, the dilepton rates from the Drell-Yan 
and charm background exceed those from the QGP. The inclusion of the correction factor 
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K^' enances the QGP dilepton yield so that it exceeds the Drell-Yan dilepton yield up to 
an invariant mass of M — 1.5 GeV, and is comparable to the charm dileptons. For the 
temperature To = 421 MeV, the QGP dilepton rate, without the correction factor 
exceeds the Drell-Yan dileptons up to an invariant mass of M = 4.5 GeV, and the charm 
dileptons up to 2 GeV. With the correction factor the thermal dilepton yield exceeds 
the DY yield up to 5.5 GeV and the charm dileptons up to 3.5 GeV, so that detection 
prospects increase considerably. 

While different temperatures have been proposed for the initial state of the QGP, the 
observation of the initial temperature is an experimental question. In this regard, higher- 
order QCD corrections, as accounted for by using the factor, enhance the possibility of 
detecting dileptons from the quark-gluon plasma for both Au+Au and S+S collisions, even 
if the initial temperature turns out to be lower than about 350 MeV. On the other hand, 
without the QCD corrections, the thermal dilepton yield from the QGP is lower than that 
from the Drell-Yan and charm decay processes and will be difficult to detect. It may be 
noted that recent calculations, approaching the first order corrections in a different manner, 
using thermal masses and finite temperature QCD effects also indicate an enhanced dilepton 
yield from the thermal plasma . 



In the plasma, the color charge of the constituents is subject to Debye screening which is 
characterized by the Debye screening length Ad, which depends on the temperature T. The 
Debye screening length at a temperature of 400 MeV is about ~ 0.2 fm. On the other 
hand, the electromagnetic annihilation into dileptons occurs within an extremely collapsed 
space zone characterized by the linear q-q distance ~ <y/y/s, which is about 0.0029 fm for 
the annihilation of a light quark- ant iquark pair at 0.5 GeV, and is much smaller than the 
Debye screening length A^. The interaction between the quark and the antiquark is not 
expected to be much affected by Debye screening in the region where annihilation occurs. 
Therefore, we expect that our correction factor for dilepton production by qq annihilation 
in the plasma will not be modified much by the addition of Debye screening corrections. 
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In the charm estimates, the open-charm production was based on a model of gluon fusion 
and qq annihilation, with an empirical K factor independent of the color state of the reacting 
partons. In this region of cc production near the threshold, it is expected that the initial- 
and final-state interactions of the participating gluons, quarks, or antiquarks are important, 
and the effective interaction depends on the color multiplet of the participating partons [IS 



A careful examination of these known effects to study open-charm production will be of 
interest to augment the investigation here and to reexamine the results of ||12|| . 

Although we have specifically focussed on the experimental situation expected at RHIC, 
the use of the correction factor is not restricted to this. In fact, the use of the can 
be extended to any experiment designed to investigate the QGP formation or existing data 
and it can be applied to account for the higher order QCD corrections to the lowest order 
thermal dilepton production in a compact way. In all situations where an intial distribution 
of the quarks and anti-quarks in the thermal plasma can be assumed, the factor can be 
used to provide estimates of the thermal dilepton production, including the higher order 
corrections and thus can serve as an useful tool in determining QGP signatures. 
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FIGURES 

FIG. 1. The differential cross section sd 2 a / 'd\frdy as a function of \fr for different rapidity y 



intervals. The data are taken from the FANL-605 experiment [21j. The solid curve are the fitting 



of the lowest order QCD calculations multiplied by the correction factor K^>f\ see text. 

FIG. 2. Dilepton production rates from the collision Au+Au at RHIC energies (\/s = 200 GeV) 
and dN/dy = 2144: (a) without the correction factor K^> and (b) with the correction factor K^ l \ 

FIG. 3. Dilepton production rates in the collision S+S at RHIC energies (y/s = 200 GeV) and 
dN/dy = 231: (a) without the correction factor and (b) with the correction factor 
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